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Abstract 

The mass of the r-lepton has been measured in the decay modes r — > 3irv T and r — > 37nr°i> T 
using a pseudomass technique. The preliminary result is 1776.71 ±0.25(stat) ±0.62(syst) MeV. The 
preliminary value of an upper limit on the relative mass difference between positive and negative 
r leptons is \(M T + - M T - ) | /M avcragc is 5.0 x 1(T 4 at 90% CL. 

PACS numbers: 
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INTRODUCTION 



Masses of quarks and leptons are fundamental parameters of the Standard Model (SM). 
In the SM high precision measurements of the mass, lifetime and leptonic branching fractions 
of the r lepton can be used to test lepton universality. Thepresent PDG value of the r 
mass ^ is dominated by the result of the BES Collaboration and has an accuracy about 
0.3 MeV. The high statistics of the Belle data allow a measurement with the same level 
of accuracy. The methods used for the r mass measurement are different for the BES and 
the Belle experiments: BES analysed the cross section for r pair production near threshold 
while in Belle the four-momenta of the visible r decay products at -^/s=10.58 GeV are used. 
This leads to different sources of systematic uncertainties. Eventually, by combining these 
two measurements we could significantly improve the accuracy of the r mass deterination. 

The analysis of individual r lepton decays allows to measure the masses of positive and 
negative r's separately and test the CPT theorem. A similar test was performed by OPAL 
at LEP [3J with the result (M T + - M T - )/M average < 3.0 x 10" 3 at 90% CL. The high statistics 
and quality of the Belle data allow us to improve this limit significantly. 

To determine the r mass we use a pseudomass technique that was first employed by the 
ARGUS 0] collaboration. This technique relies on the reconstruction of the invariant mass 
and energy of the hadronic system in hadronic r decays. 

THE PSEUDOMASS METHOD 

In a hadronic r decay (see Fig. 1) the r lepton mass M T is related to the four-momentum 
of the resulting hadronic system X by the formula 

M 2 T = M 2 X + Ml + 2E X E U - 2P X P V cos 9 ( 1 ) 

where Mx, Ex and Px are the invariant mass, energy and absolute value of the momentum 
of the hadronic system respectively; M v , E v and P v are the same quantities for the neutrino 
and 9 is the angle between the momenta of the neutrino and hadronic system. 




FIG. 1: Illustration for the variable definitions used in the Eq.(l) 

If M u = we have P v = E v = E T — Ex- This gives the following expression for the r 
mass: 

M 2 T = M 2 X + 2(E T - Ex)(E x - Px cos9). (2) 
4 



The r lepton energy E T is obtained from the energy of the electron or positron beam, 
-E'beam, in the center-of-mass (cm.) frame. All other kinematic variables listed above will 
also be evaluated in the cm. frame of the colliding beams. 

If we set the unknown value of cos# in equation (2) equal to 1, the right side of (2) will 
be smaller than the true value for Therefore, the estimator of the r mass (the so-called 
pseudomass) used in the analysis 




MC input mass 



M min (GeV) 

FIG. 2: The Monte Carlo distributions of the pseudomass M m - m for the r — > Sir^v decays and the 
input r masses equal to 1.767, 1.777 and 1.787 GeV, respectively. The fitted values of the r mass 
versus the input mass used in the MC simulation are shown together with a straight line fit. 



M min = + 2(£ beam - E X )(E X - P x ) (3) 

should be less than or equal to the r lepton mass. In the absence of initial and final 
state radiation and assuming a perfect measurement of the four-momentum of the hadronic 
system, M min is distributed below the r mass and has an edge at M T . The initial and final 
state radiation as well as the finite momentum resolution of the detector smear the shape of 
the edge for M m i n around M T . The contributions from background processes in the selected 
t + t~ sample have smooth behavior near the r mass. We can use the threshold position 
obtained from the fit to the experimental M m \ n distribution as an estimator of the r mass. 

To illustrate this method in the Belle environment we performed simulations of r decays 
into three charged pions and a neutrino with three different input r masses: the nominal 
PDG value 1.777 GeV, 1.767 GeV and 1.787 GeV0. The generated events were passed 
through the full Belle detector response simulation and reconstruction procedures. The 
resulting M min distributions for the three input r masses given above are shown in Fig. 
2. The inset in Fig. 2 shows the dependence of the fitted masses obtained from the M m j n 
distributions on the input masses used in the MC simulation. The result of the fit of this 
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dependence to the linear function / = a® + a\ x x gives = (0.1 ± 0.2) x 10 2 GeV and 
ai = 1.000 ±0.001. 

After obtaining the value of the threshold position from the fit to the experimental M m [ n 
distribution we correct it by the value, obtained from the Monte Carlo, which is equal to 
the difference between the true input value of the r mass used in the MC simulation and 
the threshold position obtained from the simulated data. 

EXPERIMENTAL PROCEDURE 

Analysis presented here is based on the data taken with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. The total integrated luminosity used in the analysis is 253 
fb" 1 . 

A detailed description of the Belle detector is given elsewhere We mention here only 
the detector components essential for the present analysis. 

Charged tracks are reconstructed from hit information in a central drift chamber (CDC) 
located in a 1.5 T solenoidal magnetic field. The z axis of the detector and the solenoid are 
along the positron beam direction, with positrons moving in the — z direction. The CDC 
measures the longitudinal and transverse momentum components (along the z axis and in the 
rip plane, respectively). Track trajectory coordinates near the collision point are provided by 
a silicon vertex detector (SVD). Photon detection and energy measurements are performed 
with a CsI(Tl) electromagnetic calorimeter (ECL). Identification of kaons is based on the 
information from the time-of-flight counters (TOF) and silica aerogel Cherenkov counters 
(ACC) . The ACC provides good separation between kaons and pions or muons at momenta 
above 1.2 GeV. The TOF system consists of a barrel of 128 plastic scintillation counters, 
and is effective in K/n separation mainly for tracks with momentum below 1.2 GeV. The 
lower energy kaons are also identified using specific ionization (dE/dx) measurements in the 
CDC. Identification of electrons is made using combined information from ECL, ACC, TOF 
and CDC 0. The magnet return yoke is instrumented to form the Kl and muon detector 
(KLM), which detects muon tracks 0] and provides trigger signals. 

The signal events were efficiently triggered by several kinds of track triggers that require 
two or more CDC tracks with combinations of TOF hits, ECL clusters or its energy sum. 
Here, we do not eliminate any events using the trigger condition information. The trigger 
conditions are complementary to each other for the detection of four-prong events in the 
present case. We can expect a high trigger efficiency (~ 95%) by combining them. 

We used only on-resonance data because the absolute beam energy calibration is known 
for this data sample better than for the off- resonance data taken 60 MeV below the T(4S) 
(y/s = 10.58 GeV ). 

t + t~ events were selected where one r lepton decays leptonically into lviv T . The other r 
lepton decays into hadronic decay modes with 3 charged pions and a neutrino or 3 charged 
pions with 1 neutral pion and a neutrino. 

The preselection of events is based on the following criteria: 

(a) Visible reconstructed energy E sum > 0.1Sy/s; 

(b) Number of well-reconstructed charged tracks greater than 2; 

(c) The sum of the z components of each good track and good photon momenta is required 
to satisfy \P Z \ < 0.5y/s. 
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Conditions (a) and (c) are calculated in the cm. frame. 
The criteria for good charged tracks are: 

• p T > 100 MeV; 

• Impact parameters Ar < 2 cm, \Az\ < 4 cm. 

Good photons are defined as ECL clusters with energy greater than 100 MeV that are 
not associated with charged tracks. The angular acceptance for photons is 17° < 6 < 150°. 
The preselection cuts suppress Bhabha, /i + \T and two-photon events. 
After the preselection the following cuts were applied. 

• Total charge equal to zero; 

• Number of leptons (muons or electrons) equal to one; 

• Number of charged pions equal to three; 

• Number of charged kaons and protons equal to zero; 

• Number of Kg's equal to zero; 

• Number of 7r°'s equal to zero for the r — > 37r^ T decay mode and equal to one for the 
r — > 37T7r°z/ r mode. 




P, = 1777.41 ± 0.25 MeV 



M min (GeV) 

FIG. 3: The pseudomass distribution M m j n for the r — ► 3^1/ decays. The points with error bars 
are data and the solid line is the result of the fit with function (4). 

The M min distribution for the r — > 37iu data is shown in Fig. 3. A fit was performed to 
the data with the function 
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FIG. 4: The pseudomass distribution M m i n for the r — > 37r decays. The points with error 
bars are data and the solid line is the result of the fit with function (4). 



F(x) = (p 3 + Pa x x) x arctan((x — p\)/p2) + Pb + Pe x z (4) 

The value of the parameter p\ obtained from the fit is p\ = 1777 .41 ± 0.25 MeV. 

The difference between the threshold position obtained from using equation (4) and the 
true value of the r mass obtained from MC is equal to Spx = 0.70±0.40 MeV. The uncertainty 
in Spi is dominantly due to limited Monte Carlo statistics (which is about 1/2 of the data) 
and the systematics of the fit procedure (choice of the fit range, shape of the threshold and 
background function). 

The distribution of M min for the 37r7r°z/ T decay mode is shown in Fig. 4 together with the 
results of the fit with the same function. 

The value of the r mass estimator for this decay mode is p\ = 1777.22 ± 0.56 MeV which 
is consistent within errors with the result from the r — > 37r ± z/ decay mode. As the statistical 
error for the "i^^v mode is significantly larger than for the 37r ± z/ one, we will concentrate 
on the former decay mode only. 



SYSTEMATIC UNCERTAINTIES 

The following contributions to the overall systematic uncertainty were considered: 

• Calibration of the tracking system. We used muons from the decay of Y(IS'), which 
are decay products of T(2S, 3S) to T(lS')7r + 7r _ . The peak position of T(1S) is shifted 
from the nominal PDG value by —4.5 ± 2.3 MeV. 

The invariant mass distribution of yU + /i~ for the T(1S) peak is shown in Fig. 5. As 
the T(15) is produced almost at rest we use the relative mass shift of the visible 
peak position from the PDG value as an estimate of the accuracy of the tracking 
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FIG. 5: The invariant mass distribution for Y(1S) — * 



calibration, i.e. AM? /My = AP^/P^. The next assumption is that the relative shift 
of the momentum of the three pion system is equal to the relative shift of the single 
track momentum, i.e. AP 3n /P 3w = AP^/P^. This allows us to propagate the shift 
of P 37T to a shift of M min . In order to estimate the shift in the threshold position we 
finally take P^ n ~ P r and M min « M T , because this is the region that is sensitive to 
the t mass. According to the above arguments, the relative shift of the T(1S) mass, 
which is equal to 4.5/9460.3 gives a systematic error from tracking of 0.39 MeV. 

• Choice of the fit range and the shape of the threshold function for the r mass es- 
timation. In addition to the function arctan((a; — pxj/p-z), where p\ and P2 stand 
for the mass position and resolution respectively, we also tried the parametrizations 
(x — pi)/ \jp2 + (x — pi) 2 and 1/ (1 + exp((x — Px)/p2}). To esimate the value of the 
correction to our estimator of the r mass and its uncertainty, we used a Monte Carlo 
t + t~ sample with a statistics of approximately one half compared to the data sample 
and an input r mass of 1777.0 MeV. We performed fits to the data and MC samples 
using three different fit functions mentioned above in five ranges of M min . 

From the fit to the MC distribution by Eq. (4) we get a difference between the visible 
threshold position and the true input r mass of Spi = 0.70 MeV. The variation of the 
fit ranges and threshold function gives a variation of the 8p\ within a ± 0.40 MeV 
range. 

We correct the estimator value p\ obtained from the fit by equation (4) to the data by 
0.70 MeV and take the value of 0.40 MeV as an estimate of the systematic uncertainty. 

• Uncertainty in the beam energy. 

For the estimation of this uncertainty we used the internal Belle analyses of the full 
reconstructed B decays for the energy calibration. In these analyses the reconstructed 
B meson energies were compared with the beam energies supplied by KEKB. The 
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conclusion from these analyses is that the beam energy is known with accuracy better 
than 1.5 MeV. This uncertainty can be translated to an uncertainty in the r mass 
from the following formula 

<r(M min ) = Ex ~ Px a(E hcam ). (5) 
ivi m - m 

Near the threshold in the M min distribution we can set E x ~ E he3im , M x ~ M T and 
M min ps M r which gives a(M) ps 0.17a(E b eam)- For a(Ebeam) = 1-5 MeV we find 
a(M T ) ps 0.26 MeV. 

To check the relation cr(M) ps 0.17a(-Ebeam), we performed a simulation of r decays 
with different -Ebeam values. We repeated the fit procedure for all MC samples assuming 
-Ebeam = M(T(AS))/2 and plotted the fit results versus the beam energy used in the 
simulations in Fig. 6. The result of the straight line fit to the plot in Fig. 6 gives 
the value of the slope of P 2 = 0.1753 ± 0.0002, which is consistent with the analytical 
calculation. 




5.28 5.2825 5.285 5.2875 5.29 5.2925 5.295 5.2975 5.3 



E be o m (GeV) 

FIG. 6: The dependence of the fit value of the threshold position on the beam energy used for the 
event simulation. The line is a result of a straight line fit to this dependence. 

• Systematic uncertainties coming from misidentified r decays are negligible, since their 
M min distributions show no significant structure in the region of the r mass. In Fig. 
7 the MC M min distributions are shown for correctly identified pions and incorrectly 
identified particles. These distributions are obtained from the same MC sample. The 
background from non-r + r~ events can also be neglected. 

Adding all these uncertainties in quadrature results in a total systematic error of 0.62 
MeV. The final result is M T = 1776.71 ± 0.25(stat) ± 0.62(syst)MeV 
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FIG. 7: The MC distribution of the pseudomass M m [ n for the decays r — > SyT^f with correctly 
identified pions. The filled histogram shows the contribution from incorrectly identified particles. 



CPT TEST 

The pseudomass method allows a separate measurement of the masses of the positively 
and negatively charged r leptons. 

A mass difference between positive and negative r leptons would result in a difference 
in the energy between the r's produced in the e + e~ collision. This in principle makes the 
assumption E T = -Ebcam invalid. The distributions of the M min for positive and negative r's 
decaying into 37rz/ are shown in Fig. 8 together with the results of the fit. 

Good agreement between the distributions for r + and r" is seen. The mass difference 
obtained from the independent fits to these distributions is M T + —M T - = —0.12 ±0.45 MeV. 

Most sources of systematic errors affect the result for positive and negative r leptons in 
the same way, so that their contributions to the mass difference cancel. However, particles 
and antiparticles interact differently with the detector material. 

To estimate a systematic shift in the measurement of particle and antiparticle momenta 
we compared the peak positions of D° — > K~tt + and D° —>■ K + tt~ , A c — > pK~n + and A c — > 
pK + 7r~, D + -> 0(lO2O)vr + and D~ -> 0(1020)71"- D s -> 0(1O2O)tt + and D s -> 0(lO2O)7r _ . 

The average relative mass shift from the decay modes listed above is about 0.8 x 10 -4 , 
which gives a systematic uncertainty in the mass difference between r + and t~ of 0.15 MeV. 

Adding the statistical and systematic errors in quadrature we obtain M T + — M T - = 
-0.12 ±0.47 MeV. 

This result can be expressed as an upper limit on the relative mass difference |§| 

\(M T + - M r -)|/M average < 5.0 x 10" 4 at 90% CL. (6) 
Without assuming CPT invariance it is no longer obvious that the charges and masses of 
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FIG. 8: The distributions of the pseudomass M m \ n for the decays r — > Sir^v separately for positive 
and negative r decays. In blue the distribution for positive, in black for negative r decays is shown. 
The solid lines are the results of the fit with function (4). 

positive and negative r decay products should be the same. Good agreement of the M min 
distributions for positive and negative r — > 37rz/ decays shows that at the present level of 
experimental accuracy CPT invariance is respected. 

RESULTS 

We have measured the mass of the r lepton from the pseudomass distributions of r decays 
into three charged pions and neutrino. The result is 

M T = 1776.71 ± 0.25(stat) ± 0.62(syst) MeV. (7) 

We obtained an independent measurement of the positive and negative r mass. The 
measured values are consistent and an upper limit on the relative mass difference is | (M T + — 
M r -)|/M average is 5.0 x 10~ 4 at 90% CL. 
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